Design, fabrication and operation of antennas for diffusive environments

ABSTRACT

A method of designing, fabricating and operating antennas is disclosed that considers the diffusive nature of the environment in which the antennas are to operate. Furthermore, the antennas can be designed, fabricated and operated so as to provide the optimal channel capacity possible given the diffusive nature of the environment in which they are to operate. The illustrative embodiment of the present invention comprises: describing an environment; describing a candidate antenna; determining a performance characteristic based on the candidate antenna with respect to the environment; and fabricating a first antenna in accordance with the candidate antenna.

REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application No. 60/125,162, filed Mar. 19, 1999, which is also incorporated by reference.

FIELD OF THE INVENTION

The present invention relates to the design, fabrication and operation of antennas in general, and, more particularly, to a technique for designing, fabricating and operating antennas that considers the diffusive nature of the environment in which the antennas are to operate.

BACKGROUND OF THE INVENTION

Although it is not usually difficult to design a workable antenna, it is notoriously difficult to design a good antenna. And while the truthfulness of this statement may be clear for the radio amateur, it is also true for the professional antenna designer who has experience, a state-of-the-art laboratory, and a modern computer with good antenna modeling software.

One of the reasons that a good antenna is difficult to design is that the elements of the antenna interact in a complex and nonlinear manner. Recently, advances in antenna modeling software have made this consideration easier. Another reason is that the objects in the environment in which the antenna operates might scatter the transmitted signal.

FIG. 1 depicts an illustrative terrestrial environment that comprises: transmitting antenna 101, receiving antenna 102, forest 111, building 112, building 113 and boat 114. As a signal is transmitted from transmitting antenna 101 to receiving antenna 102, the signal is likely to be scattered by objects in the environment that are near and between the transmitting antenna and the receiving antenna. A good antenna design considers the scattering of the transmitted signal.

In the prior art, the multipath character of the environment has not, in general, been considered in designing antennas. Rather, designers have usually made the simplifying assumption that the antennas operate in “free space.” FIG. 2 depicts a transmitting antenna and a receiving antenna in free space. When antennas are operating in free space, it is assumed that the transmitted signal radiates without scattering from the transmitting antenna to the receiving antenna. This assumption is perhaps reasonable for terrestrial microwave and satellites, but is untenable for many terrestrial applications (e.g., cities, etc.). The result is that antennas designed and fabricated to operate in free space provide poor performance when operating in diffusive environments. Therefore, the need exists for a technique for designing and fabricating antennas that considers the multipath character of the environment in which the antennas are to operate.

SUMMARY OF THE INVENTION

Some embodiments of the present invention are able to design, fabricate and operate antennas without some of the costs and disadvantages of techniques in the prior art. In particular, the illustrative embodiment of the present invention not only considers the multipath character of the environment in which the antennas will operate, but also takes advantage of the scattering to make better antennas. Furthermore, the illustrative embodiment of the present invention can design, fabricate and operate antennas that provide optimal channel capacity by taking advantage of the multipath character of the environment in which the antennas operate.

The illustrative embodiment of the present invention models the multipath character of the environment using diffusive models and uses an iterative approach to predict the performance of candidate antenna designs in that environment and to suggest improvements in the design until the predicted performance reaches an optimal or otherwise acceptable level.

The illustrative embodiment of the present invention comprises: describing an environment; describing a candidate antenna; determining a performance characteristic based on the candidate antenna with respect to the environment; and fabricating a first antenna in accordance with the candidate antenna.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 depicts an illustration of two antennas in a multipath environment.

FIG. 2 depicts an illustration of two antennas in a free space.

FIG. 3 depicts a flowchart of the illustrative embodiment of the present invention.

DETAILED DESCRIPTION

FIG. 3 depicts a flowchart of the illustrative embodiment of the present invention. First, the illustrative embodiment is described in its generalized form as it is applied to any type of antennas in any type of environment. Thereafter, the illustrative embodiment is described as it is applied to two specific examples, which are chosen to aid in an understanding of the present invention.

I. Generalized Technique

The illustrative embodiment of the present invention comprises four phases. In Phase 1 (step 301), the environment in which the antennas are to function is described. In Phase 2 (steps 302 and 303), the candidate antennas are described in terms of those parameters that if changed might affect the performance of the antennas. In Phase 3 (steps 304, 305 and 306), the performance of the candidate antennas are predicted with respect to the environment described in Phase 1. If, after Phase 3, the predicted performance is unsatisfactory, the illustrative embodiment successively iterates through Phases 2 and 3, each time varying one or more parameters of the candidate antennas, until the performance of the candidate antennas is optimal or satisfactory. In Phase 4 (step 308), the antennas are fabricated, deployed and operated in accordance with the parameters that yielded the satisfactory performance prediction.

The illustrative embodiment of the present invention predicts the performance of the antennas for a signal of interest, which by definition comprises just a single frequency defined in terms of its wavelength, λ. Antennas designed in accordance with the present invention can easily transmit and receive more than one frequency at a time, but the illustrative embodiment of the present invention only considers a signal of interest comprising one frequency at a time. It will be clear, however, to those skilled in the art how to make and use embodiments of the present invention that consider a signal of interest comprising a plurality of frequencies.

Because the illustrative embodiment of the present invention considers the nature of the environment surrounding the antennas in designing the antennas, at step 301, those aspects of the environment that might affect the propagation of the signal of interest from the transmitting antenna to the receiving antenna are described. In particular, those aspects of the environment that might affect the propagation of the signal of interest are described in terms of their properties or geometry or both.

A specific environment (e.g., Bob's Warehouse at 42nd Street and 11th Avenue, Sherwood Forest, downtown St. Louis, etc.) might be described or a nonspecific environment (e.g., a typical warehouse, a typical deciduous forest, a typical city, etc.) or a combination might be described.

The properties and geometric factors about the environment that might be described include:

Is the environment diffusive? In other words, is the mean free path of the environment much greater than the wavelength of the transmitted signal? Is all of the environment diffusive or only some portions?

If only some portions of the environment are diffusive, where are the antennas with respect to the diffusive portions? Are both the transmitting and receiving antennas deep within a diffusive portion (e.g., both within a building, one within a building and the other without, both within different buildings, etc.) or is one antenna inside a diffusive portion and the other outside the diffusive portion (e.g., the transmitting antenna is high on a tower where there is no clutter and the receiving antenna is on the ground floor of a building in Manhattan where there is lots of clutter, etc.).

Is the scattering of the transmitted signal isotropic? For example, the scattering within a building with walls at 90 degree angles is not isotropic because the scattering is not random.

Are there considerable or negligible signal losses due to absorption in the environment?

Are the signal losses due to absorption in the environment isotropic?

The way that these environmental factors can be described in a useful (i.e., quantitative) form will be described below. It will be clear to those skilled in the art what other environmental properties and geometric factors that affect the propagation of the signal of interest might be considered.

In general, there is a trade-off between considering many properties and geometric factors and ignoring the properties and geometric factors. The consideration of many properties and geometric factors of the environment will tend to:

1. increase the performance of the resulting antennas;

2. increase the computational complexity of the process for designing the antennas; and

3. decrease the interval during which the parameters chosen in accordance with the illustrative embodiment are accurate (because the environment may change over time).

Therefore, it will clear to those skilled in the art that, in general, the environment in which the antennas are to operate should probably not always be described in infinitesimal detail, but that certain simplifying assumptions should often be made. In many cases, the intentional and careful omission of some details will not affect the ability of the illustrative embodiment to design the antennas.

At step 302, the antenna elements in both the transmitting antenna and the receiving antenna are described in terms of their properties or geometric factors or both. Advantageously, the properties and geometric factors of each antenna element are described in terms of parameters that, if changed, might improve the performance of the resulting antennas.

The properties and geometric factors about the antenna elements that might be described include:

Are the antenna elements directional or omnidirectional?

What is the size of the antenna element as compared to the wavelength of the signal of interest?

What is the three-dimensional shape of the antenna element?

Does the antenna element distort the near-field signal significantly?

How much does the antenna element feed influence the signal characteristics?

The way that the properties and geometric factors of the antenna elements factors can be described in a useful (i.e., quantitative) form will be described below. It will be clear to those skilled in the art what other properties and geometric factors that affect the propagation of the signal of interest might be considered.

As in step 301, there is a trade-off between considering many properties and geometric factors and ignoring the properties and geometric factors. The consideration of many properties and geometric factors of the antennas will tend to:

1. increase the performance of the resulting antennas;

2. increase the computational complexity of the process for designing the antennas; and

3. decrease the interval during which the parameters chosen in accordance with the illustrative embodiment are accurate (because the environment may change over time).

Therefore, it will clear to those skilled in the art that, in general, the antenna elements should probably not always be described in infinitesimal detail, but that certain simplifying assumptions should often be made. In many cases, the intentional and careful omission of some details will not affect the ability of the illustrative embodiment to design the antennas.

At step 303, if either the transmitting antenna or the receiving antenna comprises a plurality of elements (i.e., is a compound antenna), the compound nature of the antennas are described in terms of their properties or geometry or both. Furthermore, the position of the antennas with respect to the environment and with respect to each other is described. Advantageously, the properties and geometry of the compound nature of each antenna are described in terms of parameters that, if changed, might improve the performance of the resulting antennas.

The properties and geometric factors about the compound nature of the antennas that might be described include:

How many antenna elements are in the transmitting antenna? How many antenna elements are in the receiving antenna? For the purposes of this specification, the number of antenna elements in the transmitting antenna is represented by n_(T), the number of antenna elements in the receiving antenna is represented by n_(R), and m=minimum(n_(T), n_(R)).

What is the geometry of the antenna elements in the transmitting antenna and in the receiving antenna? Are the antenna elements in a line? Or arranged in a two- or three-dimensional array?

Is the mutual coupling between the antenna elements to be considered or ignored?

What is the distance between the antenna elements in the transmitting antenna? What is the distance between the antenna elements in the receiving antenna? For the purposes of this specification, the distance between two antenna elements, antenna element a and antenna element b, in a single antenna is defined as r_(ab).

How are the antenna arrays pointed with respect to the environment? Up? Down? Sideways?

Is the transmitted signal power at the various transmitting antenna elements constrained or unconstrained? If it is unconstrained, the illustrative embodiment can determine the optimal distribution of power among the various transmitting antenna elements. If it is constrained, is the power evenly or unevenly distributed among the various transmitting antenna elements. For the purposes of this specification, the distribution of power among the n_(T) transmitter elements for the signal of interest is described by the transmitter power correlation matrix, M, where the trace of M equals n_(T), the matrix element M_(ij)=<x_(i)(t)x_(j) ^(*)(t)), x_(i)(t) is the normalized instantaneous signal (electric field) transmitted by transmitter element i, for i=1 to n_(T), x_(j) ^(*)(t) is the complex conjugate of the normalized instantaneous signal (electric field) transmitted by transmitter element j, for j=1 to n_(T), and <X> is the time average of X at the frequency of interest. If the transmitter power correlation matrix, M, is unconstrained, the illustrative embodiment of the present invention will compute the optimal transmitter power correlation matrix, M, in step 306.

What is the total average power at the receiving antenna from all of the transmitter elements? For the purposes of this specification, the total average power at the receiving antenna from all of the transmitting antenna elements is defined as S.

What is the noise at each receiving antenna element? In the illustrative embodiment of the present invention, the noise at each receiving antenna element is assumed to be Gaussian, independent of and identically distributed with respect to the noise at the other receiving antenna elements and its average power is assumed to be N. It will be clear to those skilled in the art how to make and use embodiments of the present invention in which the noise is not independent or identically distributed.

What is the signal to noise ratio at each receiving antenna element? For the purposes of this specification the signal to noise ratio at each receiving antenna element defined as $\rho = {\frac{S}{N}.}$

As in steps 301 and 302, the consideration of many properties and geometric factors of the compound nature of the antennas will tend to:

1. increase the performance of the resulting antennas;

2. increase the computational complexity of the process for designing the antennas; and

3. decrease the interval during which the parameters chosen in accordance with the illustrative embodiment are accurate (because the environment may change over time).

Therefore, it will clear to those skilled in the art that, in general, the compound nature of the antennas should probably not always be described in infinitesimal detail, but that certain simplifying assumptions should often be made. In many cases, the intentional and careful omission of some details will not affect the ability of the illustrative embodiment to design the antennas.

At step 304, the process of predicting the performance of the antennas described in steps 302 and 303 begins. First, certain statistical properties (e.g., the covariance, etc.) of the signal between the transmitting antenna and the receiving antenna are determined. The signal is described by G, an n_(T) by n_(R) matrix in which each matrix element G_(iα) is the signal at receiving antenna element α transmitted from transmitting antenna element i. G is a random matrix with 0 average and its covariance K is defined in terms of G as K_(ijαβ)={overscore (G_(iα)G_(jβ) ^(*))} where

G_(jβ) ^(*) is the complex conjugate of the matrix element G_(jβ) of G;

The overbar indicates an average over the multipath environment (disorder);

K is a four-dimensional matrix of size n_(T) by n_(T) by n_(R) by n_(R), comprising of elements K_(ijαβ);

for α=1 to n_(R) and β=1 to n_(R); and

for i=1 to n_(T) and j=1 to n_(T).

It will be clear to those skilled in the art how to make and use embodiments of the present invention where G has a non-zero average. In general: $\begin{matrix} {K_{{ij}\quad {\alpha\beta}} = {\int{{\hat{k}}{\int{{{\hat{k}}^{\prime}}{T_{ij}\left( \hat{k} \right)}\frac{S\left( {\hat{k},{\hat{k}}^{\prime}} \right)}{n_{T}}{R_{\alpha\beta}\left( {\hat{k}}^{\prime} \right)}}}}}} & (1) \end{matrix}$

where:

s({circumflex over (k)},{circumflex over (k)}′) is the power received at the receiving antenna from direction {circumflex over (k)}′ that is transmitted by the transmitting antenna in the direction {circumflex over (k)}, and S=∫d{circumflex over (k)}∫d{circumflex over (k)}S({circumflex over (k)},{circumflex over (k)}′);

T({circumflex over (k)}) is an n_(T) by n_(T) matrix, called the transmitter correlation matrix, in which the matrix element T_(ij) ({circumflex over (k)} is the correlation of the signal transmitted from transmitting antenna element i in the direction {circumflex over (k)} with respect to the signal transmitted from transmitting antenna element j in the same direction, and is defined as ${T_{ij}\left( \hat{k} \right)} = {\sum\limits_{\hat{e}}\quad {{w^{T}\left( {\hat{k},\hat{e}} \right)}{\chi_{i}^{T}\left( {\hat{k},\hat{e}} \right)}{\chi_{j}^{T^{*}}\left( {\hat{k},\hat{e}} \right)}}}$

where:

χ_(i) ^(T)({circumflex over (k)},ê) is the response of transmitting antenna element i to an outgoing plane wave with direction {circumflex over (k)} and polarization ê;

χ_(j) ^(T)({circumflex over (k)},ê) is the complex conjugate of the response of transmitting antenna element j to an outgoing plane wave with direction {circumflex over (k)} and polarization ê,

w^(T)({circumflex over (k)},ê) is a weight function that gives the incident power leaving in direction {circumflex over (k)} and polarization ê (where the overall scale of w^(T)({circumflex over (k)},ê) is chosen so that the trace of matrix T equals n_(T)); and $\sum\limits_{\hat{e}}$

is the normalized sum over all polarizations;

R({circumflex over (k)}′) is an n_(R) by n_(R) matrix, called the receiver correlation matrix, in which the matrix element R_(αβ)({circumflex over (k)}′) is the correlation of the signal received from receiving antenna element α from the direction {circumflex over (k)}′ with respect to the signal received from receiving antenna element β from the same direction, and is defined as ${R_{\alpha\beta}\left( \hat{k} \right)} = {\sum\limits_{\hat{e}}\quad {{w^{R}\left( {\hat{k},\hat{e}} \right)}{\chi_{a}^{R}\left( {\hat{k},\hat{e}} \right)}{\chi_{\beta}^{R^{*}}\left( {\hat{k},\hat{e}} \right)}}}$

where:

χ_(α) ^(R)({circumflex over (k)},ê) is the response of receiving antenna element α to an incoming plane wave with direction {circumflex over (k)} and polarization ê,

χ_(β) ^(R*)({circumflex over (k)},ê) is the complex conjugate of the response of receiving antenna element β to an incoming plane wave with direction {circumflex over (k)} and polarization ê,

w^(R)({circumflex over (k)},ê) is a weight function that gives the incident power arriving from direction {circumflex over (k)} and polarization ê (where the overall scale of w^(R)({circumflex over (k)},ê) is chosen so that the trace of matrix R equals n_(R)); and $\sum\limits_{\hat{e}}$

is the normalized sum over all polarizations;

∫d{circumflex over (k)} is the integral over all directions {circumflex over (k)}, normalized such that ∫d{circumflex over (k)}=1; and

∫d{circumflex over (k)}′ is the integral over all directions {circumflex over (k)}′, normalized such that ∫d{circumflex over (k)}′=1.

For isotropically diffusive environments, equation (1) becomes: $\begin{matrix} {K_{{ij}\quad {\alpha\beta}} = {R_{\alpha \quad \beta}\frac{S}{n_{T}}T_{ij}}} & (2) \end{matrix}$

where the matrix element R_(αβ) is:

R _(αβ) =∫d{circumflex over (k)}R _(αβ)({circumflex over (k)})  (3)

where ∫d{circumflex over (k)} is the integral over all directions, normalized such that ∫d{circumflex over (k)}=1. The matrix element T_(ij) is:

T _(ij) =∫d{circumflex over (k)}T _(ij)({circumflex over (k)})  (4)

where ∫d{circumflex over (k)} is the integral over all directions, normalized such that ∫d{circumflex over (k)}=1.

At the end of step 304, the covariance, K, or equivalently T({circumflex over (k)}), R({circumflex over (k)}) and S({circumflex over (k)},{circumflex over (k)}′), have advantageously been determined.

For ease of computation, it may be convenient to use an alternate basis, such as spherical harmonics, in place of direction {circumflex over (k)} and polarization ê. It will be clear to those skilled in the art that other choices of bases can be made without departing from the present invention

At step 305, a performance characteristic for the signal of interest between the receiving antenna and the transmitting antenna is determined, and, if the transmitter power correlation matrix, M, is constrained, at step 306, the value of M that optimizes the performance characteristic is determined. Advantageously, the performance characteristic is measured in terms of the channel capacity, C. In general, C is found from K— or equivalently T({circumflex over (k)}), R({circumflex over (k)}) and S({circumflex over (k)},{circumflex over (k)}′)—M, and the average of G. Here we will assume that G has 0 average. It will be clear to those skilled in the art how to make and use embodiments of the present invention where G has a non-zero average. It will be clear to those skilled in the art how to determine other performance characteristics for the signal of interest between the receiving antenna and the transmitting antenna is determined.

In the illustrative embodiment, we chose G to be known to the receiving antenna but not to the transmitting antenna. This is accomplished, for example, by having the transmitting antenna sending training sequences, periodically or sporadically, to the receiving antenna. It will be clear to those skilled in the art how to generalize this to other cases.

To reduce the computational complexity of the illustrative embodiment, there are advantageously two methods that can be used to compute the channel capacity, C, and the transmitter power correlation matrix, M. The first method is advantageously used when m is large and its accuracy is asymptotically correct as m→∞. When m is large, certain simplifying assumptions can be made that do not greatly affect the determined value of C.

The second method is advantageously used when m is small and uses Monte Carlo simulation, which is well known to those skilled in the art. In accordance with the second method, the accuracy of the determined value of C increases asymptotically with the number of Monte Carlo trials applied. The first method and the second method shall each be described in turn.

1. The First Method For Computing C and M

In general, for large m, the channel capacity, C, is found from: $\begin{matrix} {C = {{\frac{1}{\ln \quad 2}\left( {{Tr}\left\{ {\ln \left\lbrack {I_{n_{T}} + {\frac{1}{N}{\int{{\hat{k}}{\int{{{\hat{k}}^{\prime}}{T\left( \hat{k} \right)}{{MS}\left( {\hat{k},{\hat{k}}^{\prime}} \right)}{Q\left( {\hat{k}}^{\prime} \right)}}}}}}} \right\rbrack} \right\}} \right)} + {\frac{1}{\ln \quad 2}\left( {{{Tr}\left\{ {\ln \left\lbrack {I_{n_{R}} + {\frac{1}{N}{\int{{\hat{k}}{P\left( \hat{k} \right)}{R\left( \hat{k} \right)}}}}} \right\rbrack} \right\}} - {m{\int{{\hat{k}}{Q\left( \hat{k} \right)}{P\left( \hat{k} \right)}}}}} \right)}}} & (5) \end{matrix}$

where:

M is the transmitter power correlation matrix as defined above; and

Q({circumflex over (k)}) and P({circumflex over (k)}) are scalars that can be found from: $\begin{matrix} {{P\left( \hat{k} \right)} = {\frac{1}{m}{Tr}\left\{ {\frac{1}{N}{\int{{{\hat{k}}^{\prime}}{S\left( {{\hat{k}}^{\prime},\hat{k}} \right)}{T\left( {\hat{k}}^{\prime} \right)}{M\left\lbrack {I_{n_{T}} + {\frac{1}{N}{\int{{{\hat{k}}^{''}}{\int{{{\hat{k}}^{\prime}}{T\left( {\hat{k}}^{''} \right)}{{MS}\left( {{\hat{k}}^{''},{\hat{k}}^{\prime}} \right)}{Q\left( {\hat{k}}^{\prime} \right)}}}}}}} \right\rbrack}^{- 1}}}} \right\}}} & (6) \\ {{Q\left( \hat{k} \right)} = {\frac{1}{m}{Tr}\left\{ {{R\left( \hat{k} \right)}\left\lbrack {I_{n_{R}} + {\int{{{\hat{k}}^{\prime}}{P\left( {\hat{k}}^{\prime} \right)}{R\left( {\hat{k}}^{\prime} \right)}}}} \right\rbrack}^{- 1} \right.}} & (7) \end{matrix}$

For isotropically diffusive environments, equations (5), (6) and (7) become greatly simplified. In that case, the channel capacity, C, is found from: $\begin{matrix} {C = {\frac{1}{\ln \quad 2}\left\lbrack {{- {mPQ}} + {{Tr}\quad {\ln \left( {I_{n_{T}} + {\rho \quad {QTM}}} \right)}} + {\sum\limits_{a = 1}^{n_{R}}\quad {\ln \left( {1 + {PR}_{a}} \right)}}} \right\rbrack}} & (8) \end{matrix}$

where:

R_(α) is the αth eigenvalue of matrix R; $\begin{matrix} {P = {\frac{1}{m}{Tr}\rho \quad \frac{TM}{I_{n_{T}} + {\rho \quad {QTM}}}}} & (9) \\ {Q = {\frac{1}{m}{\sum\limits_{a = 1}^{n_{R}}\quad \frac{R_{a}}{1 + {PR}_{a}}}}} & (10) \end{matrix}$

At step 306, the transmitted power correlation matrix Mis determined. If the transmitter power correlation matrix, M, is constrained to a predetermined and fixed value, then equations (9) and (10) are solved, simultaneously, and resulting values for P and Q are plugged into equation (8).

If the transmitter power correlation matrix, M, is unconstrained, then the eigenvalues of M that yield the optimal value of C can be determined by: $\begin{matrix} {M_{i} = {\left( {\frac{1}{\Lambda} - \frac{1}{\rho \quad {QT}_{i}}} \right) \cdot {\Theta \left( {\frac{1}{\Lambda} - \frac{1}{\rho \quad {QT}_{i}}} \right)}}} & (11) \end{matrix}$

where

T_(i) is the ith eigenvalue of matrix T;

Θ(x)=0 if x<0 and Θ(x)=x if x≧0; and

Λ is determined together with P and Q from equations (9), (10) and equation (11) $\begin{matrix} {{\sum\limits_{i = 1}^{n_{T}}\quad M_{i}} = n_{T}} & (12) \end{matrix}$

In the case of unconstrained M, equation (9) simplifies to become: $\begin{matrix} {{PQ} = {\frac{n_{T}}{m}\Lambda}} & \text{(9a)} \end{matrix}$

M can then be found from the eigenvalues of M and the unitary matrix V of the eigenvectors of the matrix T, which is defined by: $\begin{matrix} {\begin{bmatrix} T_{1} & 0 & 0 \\ 0 & ⋰ & 0 \\ 0 & 0 & T_{n_{r}} \end{bmatrix} = {V^{+} \cdot T \cdot V}} & (13) \end{matrix}$

where V⁺ is the Hermitian conjugate of V; $\begin{matrix} {M = {{V\begin{bmatrix} M_{1} & 0 & 0 \\ 0 & ⋰ & 0 \\ 0 & 0 & M_{n_{T}} \end{bmatrix}}V^{+}}} & (14) \end{matrix}$

2. The Second Method For Computing C and M

At step 305, for small m, the channel capacity, C, is found from: $\begin{matrix} {C = {{Tr}\left\{ {\log_{2}\left( {I_{n_{T}} + {\frac{1}{N}G^{+}{GM}}} \right)} \right\}}} & (15) \end{matrix}$

where G is a Gaussian random matrix with 0 average and its covariance K is defined in terms of G as K_(ijαβ)={overscore (G_(iα)G_(jβ) ^(*))} The value of C is found by generating many random values for G, in accordance with well-known Monte Carlo techniques. It will be clear to those skilled in the art how to make and use embodiments of the present invention where G has a non-zero average.

At step 306, the transmitted power correlation matrix Mis determined as part of step 305 by varying the values for M until a satisfactory or optimal value of C is found.

At step 307, the value of C from step 305 is correlated with the parameters defined in steps 302 and 303, and the decision is made whether the value of C is satisfactory. If C is satisfactory, then control proceeds to step 308; otherwise steps 302 through 306 are iteratively repeated until an optimal or satisfactory value for C is found.

At step 308, the antennas are fabricated in accordance with the parameters defined in steps 302 and 303 that correspond to the optimal or satisfactory value for C. It will be clear to those skilled in the art how to fabricate the antennas in accordance with the parameters defined in steps 302 and 303.

At step 309, the antennas are operated in accordance with the transmitter power correlation matrix, M,computed above. It will be clear to those skilled in the art how to operate the antennas in accordance with the transmitter power correlation matrix, M.

II. The First Example

The first example is a very simple and idealized example that involves the design and fabrication of two array antennas that are both within a uniformly and isotropically diffusive environment with a mean free path that is much larger than the wavelength of the signal of interest. The wavelength of the signal of interest is λ=15 cm. Both the receiving antenna and the transmitting antenna comprise 100 (i.e., n_(T)=n_(R)=100) antennas in a line, which transmit and receive omnidirectionally. The signal to noise ratio, ρ, is 100. Furthermore, the environment is assumed to have no losses due to absorption, no parasitic affects of the antennas are considered, and there is no mutual coupling between the antenna elements. The distance between the individual antenna elements in both the transmitting antenna and the receiving antenna is represented by a and is the only parameter of the antennas that has been left to be determined by the illustrative embodiment.

Because we have assumed that the antenna elements are point antennas (i.e., the antennas are affected by all polarizations of the electric field at one point), equations (3) and (4) simplify to: $\begin{matrix} {R_{a\quad \beta} = \frac{\sin \left( {\frac{2\pi}{\lambda}r_{a\quad \beta}} \right)}{\left( {\frac{2\pi}{\lambda}r_{a\quad \beta}} \right)}} & (17) \\ {T_{ij} = \frac{\sin \left( {\frac{2\pi}{\lambda}r_{ij}} \right)}{\left( {\frac{2\pi}{\lambda}r_{ij}} \right)}} & (18) \end{matrix}$

Thereafter, the eigenvalues of R and T can be computed in well-known fashion, and C, P, Q and Λ can be computed using equations (10), (11), (12) and (13) to yield (as a function of α):

TABLE 1 Values of a as a Function of C a C (cm) (bits/sec/Hz) 4.5 415.49 5.5 466.09 6.5 510.68 7.5 548.26 8.5 540.73 9.5 536.40

From Table 1, it can be seen that the greatest value of C occurs when α=7.5 cm, and, therefore, the antennas described above should be fabricated with α=7.5 cm. Furthermore, it should noted that the capacity decreases more rapidly when α becomes less than 7.5 cm as opposed to the decrease of C when α is greater than 7.5 cm. This effect becomes more pronounced when the number of antenna elements increases.

III. The Second Example

The second example is less idealized than the first and is chosen to demonstrate another facet of the illustrative embodiment. The second example involves the design of two array antennas that each comprise 50 (i.e., n_(T)=n_(R)=50) point antennas in a line, which transmit and receive omnidirectionally. One portion of the environment is uniformly and isotropically diffusive but the other is free space and the boundary between the free space and the diffusive portion is the z=0 plane. The transmitting antenna is well above the z=0 plane and the receiving antenna is well below the z=0 plane. The wavelength of the signal of interest is λ=15 cm, and the diffusive portion of the environment has a mean free path that is much larger than λ. The signal to noise ratio, ρ, is 100. Furthermore, the environment is assumed to have no losses due to absorption, no parasitic affects of the antennas are considered, and there is no mutual coupling between the antenna elements. The distance between the transmitter antenna elements is 9 cm. and the distance between the receiving antenna elements is 7.5 cm. The transmitter correlation matrix, M, is constrained and the transmitting power is evenly distributed among all of the transmitter antenna elements. The only parameter that can be varied in this example is the angle, θ, between the line of transmitting antennas and the plane z=0.

Equations (2) and (3) & (4) are used to compute the covariance. In this example, w^(T)({circumflex over (k)},ê) does not depend on ê because we have point antennas and signal is transmitted with both polarizations. It equals

w ^(T)({circumflex over (k)},ê)=−2{circumflex over (k)}·{circumflex over (z)}Θ(−{circumflex over (k)}·{circumflex over (z)})

where {circumflex over (z)} is the unit vector pointing in upward z-direction

w^(R)({circumflex over (k)},ê) does not depend on ê because we have point antennas and signal is received with both polarizations and does not depend on {circumflex over (k)} because the receiving antenna is deep within the diffusive environment. Therefore, w^(R)({circumflex over (k)},ê)=1 for all {circumflex over (k)} and ê. ${{\chi_{i}^{T}\quad \left( \hat{k} \right)} = {{e^{I\quad \frac{2\quad \pi}{\lambda}\quad {\hat{k} \cdot {\overset{\_}{r}}_{i}}}\quad {and}\quad \chi_{i}^{R}\quad \left( \hat{k} \right)} = e^{I\quad \frac{2\quad \pi}{\lambda}{\hat{k} \cdot {\overset{\_}{r}}_{i}}}}},$

where {circumflex over (k)} is the direction of the incoming wave {right arrow over (r)} is the position of the ith antenna element with respect to the first antenna element, and I={square root over (−1)}.

Computing C for several values of θ yields the data in Table 2.

TABLE 2 Values of a as a Function of C θ/π C (bits/sec/Hz) 0.0 274 0.1 273 0.2 266 0.3 248 0.4 221 0.5 199

From Table 2, it can be seen that the greatest value of C occurs when θ=0, and, therefore, the antennas described above should be fabricated and deployed with θ=0.

It is to be understood that the above-described embodiments are merely illustrative of the invention and that many variations may be devised by those skilled in the art without departing from the scope of the invention. It is therefore intended that such variations be included within the scope of the following claims and their equivalents. 

What is claimed is:
 1. A method comprising: describing an environment; describing a candidate antenna; determining a performance characteristic based on said candidate antenna with respect to said environment wherein said performance characteristic is based on T({circumflex over (k)}), R({circumflex over (k)}) and S({circumflex over (k)},{circumflex over (k)}′), wherein T({circumflex over (k)}) is an n_(T) by n_(T) matrix, in which the matrix element T_(ij)({circumflex over (k)}) is the correlation of a first signal transmitted from an transmitting antenna element i in direction {circumflex over (k)} with respect to said first signal transmitted from an transmitting antenna element j in direction {circumflex over (k)}; S({circumflex over (k)},{circumflex over (k)}′) is the power received at a receiving antenna from direction {circumflex over (k)}′ that is transmitted by a transmitting antenna in direction {circumflex over (k)}; R({circumflex over (k)}′) is an n_(R) by n_(R) matrix, in which the the matrix element R_(αβ)({circumflex over (k)}′) is the correlation of a second signal received from an receiving antenna element α from direction {circumflex over (k)}′ with respect to said second signal received from an receiving antenna element β from direction {circumflex over (k)}′; and fabricating a first antenna in accordance with said candidate antenna.
 2. The method of claim 1 wherein said environment is at least partially diffusive.
 3. The method of claim 1 wherein said first antenna is a transmitting antenna and further comprising determining a transmitting power correlation matrix, M, for said first antenna.
 4. The method of claim 3 further comprising operating said first antenna in accordance with said transmitting power correlation matrix, M.
 5. The method of claim 1 wherein said performance characteristic is based on the channel capacity, C.
 6. The method of claim 1 wherein said performance characteristic is based on a covariance, K.
 7. A method comprising: describing an environment; describing a candidate transmitting antenna and a candidate receiving antenna; and determining a performance characteristic based on said candidate transmitting antenna and said candidate receiving antenna with respect to said environment wherein said performance characteristic is based on T({circumflex over (k)}), R({circumflex over (k)}) and S({circumflex over (k)},{circumflex over (k)}′), wherein T({circumflex over (k)}) is an n_(T) by n_(T) matrix, in which the matrix element T_(ij)({circumflex over (k)}) is the correlation of a first signal transmitted from an transmitting antenna element i in direction {circumflex over (k)} with respect to said first signal transmitted from an transmitting antenna element j in direction {circumflex over (k)}; S({circumflex over (k)},{circumflex over (k)}′) is the power received at a receiving antenna from direction {circumflex over (k)}′ that is transmitted by a transmitting antenna in direction {circumflex over (k)}; R({circumflex over (k)}′) is an n_(R) by n_(R) matrix, in which the the matrix element R_(αβ)({circumflex over (k)}′) is the correlation of a second signal received from an receiving antenna element α from direction {circumflex over (k)}′ with respect to said second signal received from an receiving antenna element β from direction {circumflex over (k)}′.
 8. The method of claim 7 further comprising: fabricating a first antenna in accordance with said candidate transmitting antenna; and fabricating a second antenna in accordance with said candidate receiving antenna.
 9. The method of claim 7 wherein said environment is at least partially diffusive.
 10. The method of claim 7 further comprising determining a transmitting power correlation matrix, M, for said candidate transmitting antenna.
 11. The method of claim 10 further comprising operating said first antenna in accordance with said transmitting power correlation matrix, M.
 12. The method of claim 7 wherein said performance characteristic is based on the channel capacity, C.
 13. A method comprising: describing a candidate transmitting antenna and a candidate receiving antenna; determining a performance characteristic based on an environment, said candidate transmitting antenna and said candidate receiving antenna wherein said performance characteristic is based on T({circumflex over (k)}), R({circumflex over (k)}) and S({circumflex over (k)},{circumflex over (k)}′), wherein T({circumflex over (k)}) is an n_(T) by n_(T) matrix, in which the matrix element T_(ij)({circumflex over (k)}) is the correlation of a first signal transmitted from an transmitting antenna element i in direction {circumflex over (k)} with respect to said first signal transmitted from an transmitting antenna element j in direction {circumflex over (k)}; S({circumflex over (k)},{circumflex over (k)}′) is the power received at a receiving antenna from direction {circumflex over (k)}′ that is transmitted by a transmitting antenna in direction {circumflex over (k)}; R({circumflex over (k)}′) is an n_(R) by n_(R) matrix, in which the the matrix element R_(αβ)({circumflex over (k)}′) is the correlation of a second signal received from an receiving antenna element α from direction {circumflex over (k)}′ with respect to said second signal received from an receiving antenna element β from direction {circumflex over (k)}′; and fabricating a first antenna in accordance with said candidate transmitting antenna.
 14. The method of claim 13 wherein said environment is at least partially diffusive.
 15. The method of claim 13 further wherein the steps of describing and determining are performed iteratively until said performance characteristic is satisfactory.
 16. The method of claim 13 further comprising determining a transmitting power correlation matrix, M, for said candidate transmitting antenna.
 17. The method of claim 16 further comprising operating said first antenna in accordance with said transmitting power correlation matrix, M.
 18. The method of claim 13 wherein said performance characteristic is based on a statistical property of a signal between said candidate transmitting antenna and said candidate receiving antenna. 